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ABSTRACT 

We report on the Swift monitoring of the candidate supergiant fast X-ray transient (SFXT) 
IGR J 16418 — 4532, for which both orbital and spin periods are known (~3.7 d and ~1250s, 
respectively). Our observations, for a total of ~43 ks, span over three orbital periods and 
represent the most intense and complete sampling of the light curve of this source with a 
sensitive X-ray instrument. With this unique set of observations, we can address the nature 
of this transient. By applying the clumpy wind model for blue supergiants to the observed 
X-ray light curve, and assuming a circular orbit, the X-ray emission from this source can be 
explained in terms of the accretion from a spherically symmetric clumpy wind, composed of 
clumps with different masses, ranging from ~5 x 10 16 to 10 21 g. Our data suggest, based on 
the X-ray behaviour, that this is an intermediate SFXT. 
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1 INTRODUCTION 

The X-ray transient IGR J 16418 — 4532 was discovered by INTE- 
GRAL during observations of the black hole X-ray transient 4U 
1630—47 on 2003 February 1-5 (Tomsick et al. 2004), with a 20- 
40 keV flux of 3 x 10~ u erg cnT 2 s -1 . Sguera et al. (2006) reported 
on further INTEGRAL observations, taken on 2004 February 26, 
during which fast (~1 h) X-ray outbursts were observed that peaked 
at ~80 mCrab (20-30 keV). This behaviour led Sguera et al. (2006) 
to propose that IGR J 16418 — 4532 is a candidate supergiant fast 
X-ray transient (SFXT). SFXTs (e.g. Sguera et al. 2005) are high- 
mass X-ray binaries (HMXBs) associated with OB supergiant stars 
via optical spectroscopy that show an X-ray dynamic range of three 
to five orders of magnitude. Indeed, they are characterized by bright 
(peak luminosities of 10 36 — 10 37 erg s _1 ), short (a few hours, as ob- 
served by INTEGRAL ; Sguera et al. 2005; Negueruela et al. 2006) 
X-ray outbursts significantly shorter than those of typical Be/X-ray 
binaries, and a quiescent luminosity of ~10 32 erg s _1 (e.g. in’t Zand 
2005; Bozzo et al. 2010). 

In their systematic analysis of the INTEGRAL observations from 
2003 to 2009, Ducci, Sidoli & Paizis (2010) detected 23 outbursts, 
for an activity duty cycle of ~1 per cent (one of the highest among 
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the 14 SFXTs and SFXT candidates they examined) and fluxes 
ranging between 1.3 x 10~ 10 and 4.8 x 10~ 10 erg cnT 2 s _1 (18— 
lOOkeV). 

XMM-Newton observations (Walter et al. 2006) showed a heav- 
ily absorbed [(V H = (1.0 ± 0.1) x 10 23 cm~ 2 ] pulsar (P spin = 
1246 ± 100s) with an unabsorbed flux of 1.3 x 10~ 10 erg cm -2 s^ 1 
(2-100 keV) in a fit performed in conjunction with the INTEGRAL 
data. 

Chaty et al. (2008) proposed, based on the XMM-Newton posi- 
tion, 2MASS J16415078— 4532253 as the best near-infrared coun- 
terpart candidate; they also fitted the spectral energy distribution 
(SED), obtaining a temperature of 32 800 K for the massive com- 
panion, consistent with an OB spectral classification. SED fitting by 
Rahoui et al. (2008) yielded A v = 14.5 mag, and an 08.5 spectral 
type; depending on the luminosity class they derive a distance of 
4.9 (V), 8.3 (III) or 13kpc (I). 

A periodicity of ~3.75 d was discovered by Corbet et al. (2006) 
in the SwiftfBuvst Alert Telescope (BAT; Barthelmy et al. 2005) 
(3.753 ± 0.004 d. MJD 53360-53630) and RXTE/Ml Sky Monitor 
(ASM) (3.7389 ± 0.0004 d, MJD 5009 1-53810) light curves, which 
was interpreted as the orbital period of a binary system. The hint of 
a total eclipse at phase 0.5 observed in the BAT data is consistent 
with a supergiant companion. Levine et al. (201 1) refined the orbital 
period value to 3.73886 ± 0.00028 d, by using the RXTEI ASM data 
up to MJD 55224. The observed values of the orbital and pulse 
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periods indeed place IGR J 1 64 1 8 — 4532 in a region of the Corbet 
diagram (Corbet 1986) where wind-accreting HMXBs generally lie. 

In this paper, we analyse all the Swift data collected on IGR 
J 164 18 — 4532. The first set of observations comprises two 2ks tar- 
get of opportunity (ToO) observations performed on 201 1 February 
1 8 and March 4, following a Monitor of the All Sky X-ray Image 
(MAXI) alert of a possible detection of IGR J1 64 18— 4532 on 2011 
February 17 [preliminary results were described in Romano et al. 
(2011b)]. The second set is drawn from a monitoring campaign in 
2011 July that spans over three orbital periods. This Swift monitor- 
ing represents the most intense and complete sampling along the 
orbital period of the light curve of this source with a sensitive X-ray 
instrument. 

2 OBSERVATIONS AND DATA REDUCTION 

Table 1 reports the log of the Swift/X - ray Telescope (XRT; Burrows 
et al. 2005) observations used in this paper. The Swift observations of 
IGR J 1 64 18-4532 during 20 1 1 February and March (0003 1 92900 1 
and 0003 1 929002) were obtained (Romano et al. 20 1 1 b) as two ToO 
observations following a MAXI alert (5609524433, later retracted, 
as due to GX 340+0) of a possible detection of IGR J 1 64 1 8 — 4532 
on 2011 February 17. The data collected in 2011 July were ob- 
tained as a ToO monitoring programme of nine scheduled daily 
observations, each 5 ks long, starting on 201 1 July 13. This strategy 
was devised in order to cover the light curve in different orbital 
phases while maintaining the observing time per day reasonably 
short, thus not hindering gamma-ray burst (GRB) observations. 
However, most observations were cut short by the occurrence of 
GRBs and the discovery of a new soft gamma-ray repeater (starting 
from obsID 00031929005, see Table 1). Observation 00032037001 
was performed to determine the best offset from the intended tar- 
get in order to alleviate the problem of single-reflection rings from 
the nearby source GX 340+0 (these rings marginally affect the 
source while in the faint states, so that the analysis of the earlier 
segments was still feasible, with an appropriate choice of extrac- 
tion regions). A new pointing (hence a new obsID) was adopted 


Table 1. Summary of the Swift/XRT observations. 


ObsID 

Start time (ut) 

End time (ut) 

Exp. 

(s) 

(f) 0 

0003 1929001 6 

2011-02-18 

01:59:43 

2011-02-18 

03:43:26 

1951 

0.27 

00031929002 

2011-03-04 

00:03:19 

2011-03-04 

05:01:57 

2001 

0.01 

00043 105 00 l c 

2011-06-27 

05:14:54 

2011-06-27 

05:23:56 

537 

0.80 

00031929003 

2011-07-13 

09:46:55 

2011-07-13 

14:43:57 

4756 

0.16 

00031929004 

2011-07-14 

12:39:39 

2011-07-14 

19:31:56 

4832 

0.47 

00031929005 rf 

2011-07-15 

06:42:44 

2011-07-15 

23:00:58 

273 

0.72 

00031929005 

2011-07-15 

06:42:50 

2011-07-15 

23:01:57 

1078 

0.72 

00031929006 rf 

2011-07-16 

10:15:02 

2011-07-16 

23:05:36 

458 

0.01 

00031929006 

2011-07-16 

10:15:05 

2011-07-16 

23:06:58 

379 

0.01 

00032037001 

2011-07-18 

21:28:28 

2011-07-19 

00:51:31 

1947 

0.62 

00031929007 

2011-07-19 

01:52:41 

2011-07-19 

23:19:57 

2673 

0.77 

00032037002 

2011-07-20 

05:18:42 

2011-07-20 

15:22:58 

3483 

0.01 

00032037003 

2011-07-21 

10:06:18 

2011-07-21 

15:27:54 

3761 

0.31 

00032037004 

2011-07-22 

15:19:43 

2011-07-22 

18:44:57 

2812 

0.62 

00032037005 

2011-07-23 

11:48:19 

2011-07-23 

15:23:56 

4127 

0.85 

00032037006 

2011-07-29 

12:13:30 

2011-07-29 

23:34:57 

4650 

0.50 

00032037007 

2011-07-30 

01:01:30 

2011-07-30 

09:27:56 

4287 

0.63 


“Mean phase referred to r e p 0C h =MJD 53560.2 (Levine et al. 2011). 
^Corresponding to MJD 55610.08. 

'’Serendipitous observation. 
d WT mode data. 


afterwards. In our work, we also include observation 00043105001 
that serendipitously covers the region of IGR J 1 64 18 — 4532. The 
2011 July campaign lasted 18 d divided in 12 observations for a 
total on-source exposure of ~39 ks. 

For our binary orbital period search and analysis, we retrieved 
the BAT ‘orbit-by-orbit’ light curves (averaged over Swift’s orbital 
period of ~90m, 15-50keV) covering the data range from 2005 
February 12 to 2011 July 12 (MJD range 53413-55754) from the 
BAT transient monitor (Krimm et al. 2006, 2008) page. 1 These data 
were further screened to exclude bad quality points (quality flag 1 
and 2) and referred to the Solar system barycentre (SSB) by using 
the task earth2Sun. 

The XRT data were processed with standard procedures 
(xrtpipeline vO. 12.6), filtering and screening criteria by us- 
ing ftools. Both windowed-timing (WT) and photon-counting 
(PC) events were considered. The selection of event grades 
was 0-2 and 0-12, for WT and PC data, respectively 
(Burrows et al. 2005). Source events were accumulated within an 
annular/circular region (depending on whether pile-up correction 
was required or not, respectively) with an outer radius of 20 pixel 
(1 pixel ~2.36 arcsec), and background events were accumulated 
from source-free regions far from the contamination of the single- 
reflection rings from GX 340.0+0. For our timing analysis, we 
converted the event arrival times to the SSB with the task barycorr. 
Light curves were created for several values of signal-to-noise ratio 
(S/N) and number of counts per bin; all were corrected for point 
spread function (PSF) losses, vignetting and background. For our 
spectral analysis, we extracted events in the same regions as those 
adopted for the light-curve creation; ancillary response files were 
generated with xrtmkarf, to account for different extraction re- 
gions, vignetting and PSF corrections. We used the latest spectral 
redistribution matrices in CALDB (20110705). 

The UV/Optical Telescope (UVOT; Roming et al. 2005) observed 
IGR J 1641 8 — 4532 simultaneously with the XRT with the ‘filter of 
the day’ (FoD), i.e. the filter chosen for all observations to be carried 
out during a specific day in order to minimize the filter wheel usage. 
The exception is observation 00032037007, which was requested 
in the u filter, as opposed to FoD, to attempt to obtain a detection in 
that band. The data analysis was performed using the uvotimsum and 
uvotsource tasks included in the ftools. The latter task calculates 
the magnitude through aperture photometry within a circular region 
and applies specific corrections due to the detector characteristics. 

All quoted uncertainties are given at 90 per cent confidence level 
for one interesting parameter unless otherwise stated. The spectral 
indices are parametrized as F v oc v~“, where F v (erg cm~ 2 s' 1 Hz -1 ) 
is the flux density as a function of frequency v; we adopt T = a + 1 
as the photon index, N(E) oc E~ r (photon cm -2 s _1 keV~'). 

3 RESULTS 

3.1 Orbital period from BAT data 

For the orbital modulation, we considered the BAT data in the 
time range MJD 53414—55754. With a standard folding analysis of 
the SSB light curves, we measured an orbital period of 3.740 ± 
0.002 d (lcr error), which is fully consistent with the value obtained 
by Levine et al. (2011), P orb = 3.73886 ± 0.00028 d. Therefore, 
hereafter we shall adopt the more precise measure by Levine et al. 
(2011), and we shall refer to Tepoch =MJD 53560.20000 (Corbet 

1 http://swift.gsfc.nasa.gov/docs/swift/results/transients/ 
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Figure 1. BAT light curve folded at P OI b = 3.73886 d and T e poch =MJD 
53560.20000. 

et al. 2006). Fig. 1 shows the BAT light curve folded over this 
period. 

The BAT data show an eclipse centred at <p ~ 0.55. The depth 
is consistent with a total eclipse, and the eclipse full width at half- 
maximum duration is 0.17 ± 0.05 of the orbital period. The folded 
light curve is quite flat; the peak-to-trough pulsed fraction, PF pt = 
(F max - FnuJ/CFm ax + F m i„), where F max and F min are the observed 
background-subtracted count rates at the peak and at the minimum, 
is 78 ± 12 per cent (Iff error). 


3.2 X-ray position 

We used 20.4 ks of PC mode data and simultaneous Swift/UVOT 
images to obtain an astrometrically corrected position (see Goad 
et al. 2007; Evans et al. 2009) of RA(J2000) = 250.46104, 
Dec.(J2000) = —45.54092, which is equivalent to RA(J2000) = 
16 h 41 m 50 s .65, Dec. (J2000) = — 45 0 32'27'.'3, with an uncertainty of 
1.9 arcsec (90 per cent confidence level). 2 This position is 2.5 arcsec 
from the source 2MASS J16415078— 4532253 and 4.3 arcsec from 
the XMM-Newton position (Walter et al. 2006), and it allows us 
to confirm the Two Micron All Sky Survey source as the optical 
counterpart of IGR J 1 64 1 8 — 4532. 


3.3 UVOT optical counterpart 

At the XRT position of IGR J 1641 8 — 4532, no detection was 
achieved in any ultraviolet band down to a limit of w\ > 21.78, 
m2 > 21.66, w2 > 21.89mag, respectively. A marginal detection 
(2.9ff ) was obtained in the u band, at 22.28 ± 0.40 mag. These mag- 
nitudes are on the UVOT photometric system described in Poole 
et al. (2008), and are not corrected for Galactic or intrinsic ex- 
tinction. At this position, the column due to Galactic absorption is 
1 .59 x 10 22 cm~ 2 (Kalberla et al. 2005). By adopting the conversion 
into Ay (Predehl & Schmitt 1995), N H = 1.79 x 10 21 Ay cm~ 2 , this 
corresponds to (Cardelli, Clayton & Mathis 1989) A„ = 14.2 mag 
in the UVOT u band. 


2 This new UVOT-enhanced XRT position takes advantage of the new cal- 
ibration of the XRT to UVOT detector coordinates of 201 1 August (Evans 
2011); see also http://www.swift.ac.uk/reprocessed.php. 
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3.4 Spin period search 

We looked for evidence of the 1 246 ± 1 00 s spin periodicity reported 
by Walter et al. (2006) in the SSB XRT light curves. We only 
considered observations where the source was significantly detected 
(thus excluding obsIDs 00032037001 and 00032037004). 

We searched each single observation longer than one spin period 
and the total combined data set within ±300 s [corresponding to 
a ±3<r error reported by Walter et al. (2006)] by applying both 
a Fourier transform and a Z 2 analysis (a method that does not 
require the binning of the data; Buccheri et al. 1983). No statistically 
significant signal was detected either in the individual observations 
or in the combined data set. With respect to the latter, the upper 
limit on the pulsed fraction (defined as the semi-amplitude of the 
sinusoidal modulation divided by the mean count rate), computed 
according to Vaughan et al. ( 1994), is 42 per cent at the 99 per cent 
confidence level and for periods between 946 and 1546 s. 

This upper limit is only marginally compatible with the (64 ±10) 
per cent value measured by Walter et al. (2006). We note, however, 
that they were able to detect pulsations only during a fraction of 
their XMM-Newton data and in particular while the source was not 
flaring, and that our data were collected in a variety of states of 
brightness and across almost all orbital phases. So our upper limit 
and their measure may not be straightforward to compare. 

3.5 Light curves 

Fig. 2 shows the 0.3-10keV light curve of IGR J 1 64 18 — 4532 
during the 2001 February 18 observation, binned in order to achieve 
an S/N of 5. The source reached ~5 count s _1 during the initial flare 
and then decayed through a series of smaller flares. The lowest point 
is at 0.16counts _1 . When binned at an S/N of 3, the dynamical 
range, during this observation, is of the order of 50. 

Fig. 3 shows the 0.3-10keV light curve of IGR J16418 — 4532 of 
the whole 201 1 July campaign after bin-by-bin background subtrac- 
tion and corrections for pile-up, PSF losses and vignetting, binned 
at S/N = 3. The light curve starts at phase 0.16 (assuming a period of 
P orb = 3.73886 d and an initial epoch F ep och = MJD 53560.20000) 
and monitors the source state through over three full periods. Super- 
imposed on the long-term orbital modulation, which follows that 
seen in the BAT data, flaring is observed on short time-scales, as 



100 1000 10 4 
Time (s) 


Figure 2. Swift/XKT 0.3-10keV light curve of IGR J 164 18 — 4532 during 
observation 00031929001, background-subtracted, corrected for pile-up, 
PSF losses and vignetting (S/N = 5). Time is referred to the start of the 
observation (2011-02-18; 01:59:46). 
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Figure 3. Swift/X RT 0.3-10keV light curve of IGR J164 18 — 4532 during the 2011 July monitoring programme, background-subtracted and corrected for 
pile-up, PSF losses and vignetting (S/N=3). Downward-pointing arrows are 3<r upper limits. Different colours (see the electronic edition) mark different 
observations (see Table 1). The top axis reports the phase (with P or b = 3.73886 d and Tepoch = MJD 53560.20000). 




Phase 


Figure 4. Swi/r/XRT 0.3-10keV flux light curve of IGR J 1 64 18 — 4532, folded at P or b = 3.73886d and r ep och =MJD 53560.20000. The data were collected 
in 201 1 February and June (dark grey filled circles), and 201 1 July (black filled circles). Downward-pointing arrows are 3a upper limits. The open blue circle 
is a detection obtained by combining three observations (00032037001, 00032037004 and snapshots 1-3 of 00032037007). Also included are the 2004 August 
XMM data (light grey curve). 


shown in Fig. 2. This behaviour has been observed in most SFXTs 
(Sidoli et al. 2008; Romano et al. 2009b, 2011a, 2010). 

The lowest point (as observed by XRT) was collected dur- 
ing the eclipse, by combining three observations (00032037001, 
00032037004 and snapshots 3 1-3 of 00032037007; see open blue 
point in Fig. 4) at 2.5 x 10” 3 counts” 1 . This corresponds to 
an unabsorbed 0.3-10keV flux of ~6.0 x 10” 13 erg cm” 2 s” 1 . 


3 Part of an observation with a continuous viewing of the target, with typical 
time-scales of 10 3 s. 


and to a luminosity of 1.2 x 10 34 ergs”' (assuming the opti- 
cal counterpart distance of 13kpc). The lowest point during the 
campaign outside the eclipse was recorded on MJD 55763.43 
at 0.022 count s~*, corresponding to an unabsorbed 0.3-10 keV 
flux of ~4.7 x 10” 12 erg cm” 2 s' 1 , and to a luminosity of 9.5 x 
10 34 ergs”'. As the peak count rate is reached on MJD 55610.08 
at >8 counts” 1 , the observed dynamical range of this source is at 
least 370 (1400 considering the points within the eclipse). 

Fig. 4 shows the 0.3-10 keV light curve, folded using P orb = 
3.73886 d and r epoch =MJD 53560.20000. We adopted count 
rate to flux conversion factors derived from spectral fits of each 

© 201 1 The Authors, MNRAS 419, 2695-2702 
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Count rate (0.3-10 keV) 



Phase 

Figure 5. Swift/XKT HR (red open circles, 30 count bin -1 ) as a function of 
count rate (top) and phase (bottom). The black points are weighted means. 


observation (Section 3.6) and subsequent interpolation across ob- 
servations where fitting was not feasible. We also include the ~ 10 ks 
2004 August 19 XMM-NewtonfEPlC-pn observation (Walter et al. 
2006), binned at 256 s, and converted into flux by adopting a single 
conversion factor obtained from fitting the mean spectrum with an 
absorbed power-law model [T = 1.40 ± 0.05, N a = (10.5 ± 0.5) x 
10 22 cm -2 , 1 a errors; x 2 = 0.99/169 degrees of freedom (d.o.f.)]. 
The eclipse lasts A <p ~ 0.2 and is centred at <j> ~ 0.55, as observed 
in the BAT data (Fig. 1). 

Fig. 5 shows the hardness ratio (HR; defined as HIS, where H is 
the count rate in the 4-10keV band and S is the count rate in the 
0.3-4 keV band), as a function of count rate (top) and orbital phase 
(bottom), where the light curve was binned to have 30 count bin -1 
(red open circles). A fit to a constant value yields an HR value 
of 1.01 ± 0.03 and x 2 = 0.85 for 130 d.o.f., so that the HR is 
shown to be roughly constant over two orders of magnitude in 
count rate. Furthermore, the HR were also weighted across several 
phase intervals (black points); we cannot find evidence of variations 
of the HR as a function of phase. 

3.6 Spectral analysis 

We first extracted PC spectra for each segment in which a detection 
was obtained and a minimum of ~ 100 source counts were available. 

© 201 1 The Authors. MNRAS 419, 2695-2702 
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Table 2. XRT spectroscopy of IGR J16418 — 4532. 


ObsID 

Wh 

r 

Flux“ 

L b 

X 2 /d.o.f. 


(10 22 cm -2 ) 




C-stat (d.o.f.) c 

00031929001 

6 77+ 2 - 32 
o. / /_j 78 

1 "20+0. 61 
^ -0.54 

2 1 _36+1;^ 

43 

0.532/18 

00031929002 

4-i7 IS 

i 04+0.51 

1 -U — 0.45 

11. 5+}; 39 

23 

0.845/17 

00043105001 

5-28 + J 2 *l 

0.87 

o qc+1.16 

0.83 

8 

70.05(76) 

00031929003 

5 20 +1 ‘°° 
0.82 

] oc+0.27 
I.JJ_0 24 

00 
k> 
-0 

1 + 
o o 

Lg On 

17 

0.935/58 

00031929004 

5 16 +4 ' 17 
AO — 2.81 

0 09 +0 - 87 
u,u -0.75 

1 05 +0 - 21 
A - UD -0.17 

2 

169.46(141) 

00031929007 

6, qo+1.24 
-0.97 

1 04 +0 - 40 
1 -U — 0.36 

7 62 +0 ‘ 72 
1 - oz — 0.69 

15 

0.647/28 

00032037002 

4 SQ+ 1 - 23 
H ' jy -1.00 

1 05 +0 ' 38 

I.UJ— 0 34 

6 32 +0 56 
0.54 

13 

0.831/31 

00032037003 

d , dd -1.96 

0 62 +0 ' 52 
— 0.55 

9 ci +0.32 
z - ja -0.31 

5 

1.114/16 

00032037005 

5 14+ 1 - 72 

1 oi +0 - 45 

A - UA -0.40 

o oo+0.40 
-0.37 

8 

0.838/27 


“Unabsorbed 2-10 keV fluxes (10 -11 erg cm -2 s -1 ). 
i ’2-10keV luminosities in units of 10 35 ergs -1 , at 13kpc. 
“Cash statistics (C-stat). 


We did not consider the WT data for our spectroscopic study as 
they were statistically poor and were heavily contaminated by the 
background from the single-reflection rings from GX 340.0+0. 
The data were rebinned with a minimum of 20 count per energy 
bin to allow x 2 fitting, with the exception of low-count-statistics 
spectra, in which Cash (Cash 1979) statistics and spectra binned 
to 1 count bin -1 were used, instead. The spectra were fitted in the 
0.3-10 keV energy range with a single absorbed power-law model. 
More complex models were not required by the data. The results 
are reported in Table 2. We note that the spectral parameters do not 
show significant variability within the large uncertainties, despite 
the large observed variations in flux throughout the campaign. 

Observation 00031929001 (shown as the dark grey filled circles 
at phase 0.27 in Fig. 4, and detailed in Fig. 2), in particular, deserves 
special attention. Although it was performed as a follow-up of a false 
trigger (later attributed to the nearby GX 340+0), it none the less 
found the source in a very active state, the brightest of the whole 
campaign, with its average 2-10keV unabsorbed flux of ~2 x 
10~ 10 erg cm -2 s -1 , corresponding to ~4 x 10 36 erg s -1 at 13 kpc). 
The peak flux was ~2 x 10 -9 erg cm -2 s -1 (~4 x 10 37 ergs -1 ). 
An absorbed power-law model yielded a high absorbing column 
Ah = (6.8^ 2 |) x 10 22 cm -2 (in excess of the Galactic one, 1.59 x 
10 22 cm' 2 ) and a photon index T = 1.4 +q'j (see Table 2). 

4 DISCUSSION 

IGR J 16418 — 4532 is a candidate SFXT, based on the hard X-ray 
behaviour (Sguera et al. 2006), for which we know both the orbital 
(Borbit ~ 3.7 d; Corbet et al. 2006) and the spin (~ 1250 s; Walter 
et al. 2006) periods. These values place IGR J 1 64 18 — 4532 in the 
region of the Corbet diagram (Corbet 1986) of the wind-accreting 
supergiant HMXBs. In this paper, we have taken advantage of the 
short orbital period which makes IGR J 1 64 1 8 — 4532 a very good 
source to monitor with the XRT, to study the X-ray properties of 
this source as a function of the orbital phase, in particular to address 
the issue of its SFXT nature. 

The Swift/BXT and XRT folded light curves ofIGRJ16418— 4532 
(Figs 1 and 4) show the presence of a dip that is consistent with 
zero intensity. When interpreted as an eclipse, it can be used to 
infer the nature of the primary of the binary system. The BAT 
light curve of IGR J 1 64 1 8 — 4532 (excluding the eclipse) does not 
clearly show a flux modulation with the orbital phase, suggesting 
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that the eccentricity is not large. This is in agreement with the 
observation that in many X-ray binary systems with nearly circular 
orbits and short orbital periods (<4 d), the eccentric orbit produced 
by the explosion of the supernova is circularized on time-scale of 
10 6 yr (consistent with the evolution time-scale of an O star), due 
to viscous tidal interactions (see e.g. Press, Smart' & Wiita 1975). 
Thus, for simplicity, we assume a circular orbit (P OI \, = 3.73886 d). 
Let us adopt the typical values of stellar masses and radii (Martins, 
Schaerer & Hillier 2005) for a primary of spectral type 08.5 and 
luminosity classes I, III or V, and for the neutron star a mass of 
M ns = 1.4 Mg and a radius of Rns = 10 km. By combining Kepler’s 
third law with the relationship given by Rappaport & Joss (1983) 
for the radius of the primary R, 

( R/a ) 2 = cos 2 i + sin 2 i sin 2 9 e , (1) 

where a is the semimajor axis of the orbit, i is the inclination and 
9 t is the eclipse semi-angle, we find that the observed A <p = 0.2 is 
inconsistent with a luminosity class III or V primary. We therefore 
conclude that indeed the primary is an 08.51 supergiant star [with 
mass, radius, luminosity, effective temperature given by Martins 
et al. (2005): M 0 b = 31.5Mq,R OB = 21.4 Rq, logL OB /LQ = 5.65, 
T eff = 32274 K], Therefore, the distance to IGR J 1 64 18 — 4532 is 
~13 kpc (Rahoui et al. 2008). 

During the eclipse of IGR J 1 64 18 — 4532, we detect a non-zero 
flux with the XRT (L x ~ 3 x lO^ergs' 1 ; see Fig. 4). Detections 
of HMXBs during the eclipse are quite common (see e.g. Clark, 
Minato & Mi 1988; Sako et al. 1999; van der Meer et al. 2005) and 
can be ascribed to X-ray photons scattered by the stellar wind of 
the OB star and discrete fluorescent and recombination lines from 
highly ionized gas around the compact object. The decrease in X- 
ray flux, when the compact object is approaching the eclipse, is 
related to the increase of the column density along the line of sight, 
through the stellar wind produced by the OB star. From observations 
of other wind-fed X-ray pulsars, such as 4U 1700—377, Vela X-l, 
OAO 1657—415 (see e.g. Sako et al. 1999; van der Meer et al. 
2005; Audley et al. 2006), the absorption column density during the 
eclipse ranges from ?»10 23 up to ^10 24 cm' 2 . 

We can compare the X-ray luminosity and absorption column 
density of IGR J16418 — 4532 during the eclipse with those calcu- 
lated by van der Meer et al. (2005) for 4U 1700—377, a wind-fed 
HMXBs with properties similar to IGR J 164 18— 4532 (orbital pe- 
riod P orb = 3.4 d, companion star of spectral type 06.5Iaf). van der 
Meer et al. (2005) observed 4U 1700—377 with XMM-Newton dur- 
ing an eclipse; they found an X-ray luminosity of ?»7 x 10 34 erg s' 1 
(in agreement with the X-ray luminosity of IGR J 1641 8— 4532 dur- 
ing the eclipse) and a column density of ^10 24 cm' 2 . 

The XRT light curve (Figs 2 and 3) shows a large dynamic range, 
up to one to two orders of magnitude within a single snapshot 
and about three orders of magnitude overall (even when the low- 
intensity state is excluded). While this dynamic range is lower than 
that observed in the classical (Chaty, Zurita Heras & Bodaghee 
2010) SFXTs, which reach four to five orders of magnitude in dy- 
namic range (e.g. IGR J17544— 2619), it is indeed typical of inter- 
mediate SFXTs (e.g. IGR J1 8483-03 11). The 2011 February flare 
(Fig. 2 and dark grey filled circles at phase 0.27 in Fig. 4), in par- 
ticular, reached ~2 x 10' 9 erg cm' 2 s' 1 (2-10 keV, unabsorbed), 
corresponding to a luminosity of ~4 x 10 37 ergs' 1 , which is typ- 
ical of SFXT outbursts as we have observed them with Swift/XKT 
(e.g. Romano et al. 2008, 2009a; Sidoli et al. 2009b, a). 

On the other hand, while the typical SFXT light curve is character- 
ized by a faint flaring activity with luminosity of 10 33 — 10 34 erg s' 1 
(Sidoli et al. 2008; Romano et al. 2009b, 2011a), sporadically inter- 


rupted by flares with duration of a few hours and peak luminosity 
of lO^-lCPergs' 1 , IGR J 16418 — 4532 exhibits a flaring activity 
ranging from ~10 35 up to ~4 x 10 37 erg s' 1 [also see Walter et al. 
(2006) and the light grey curve in Fig. 4], 

To reproduce the observed X-ray variability, we applied the 
clumpy wind model of Ducci et al. (2009), which we briefly sum- 
marize here. This model assumes that the wind of OB supergiants 
is inhomogeneous, composed of dense clumps surrounded by a 
hotter and homogeneous wind. Clumps follow a power-law mass 
distribution: 

/M cl V f 

P(M C1 ) = /:(-] , (2) 

where k — (1 — ft) Mf* — M a '' f ) is the normalization 
constant, M c \ is the mass of a clump, and M a and M* define the 
extrema of the mass range. Assuming a spherical geometry for 
the clumps, the initial clump dimension distribution is given by 
Nm c i oc R^ clump s' 1 , where R c i is the radius of the clump. The 
clump mass-loss rate M c \ is related to the total mass-loss rate M tot 
by means of / = M cl /M tot , the fraction of wind mass in the form of 
clumps. Clumps and the homogeneous wind follow the j3 -velocity 
law v(r) = Uoo (l — 0.9983^^)^ (Castor, Abbott & Klein 1975), 
where /J is a constant in the range ~0.5-1.5, Rob is the radius of 
the supergiant and Voo is the terminal wind velocity. 

IGR J 1 64 18 — 4532 was intensely monitored by Swift/XKT for 
two orbital cycles (MJD 55755-55766). Each snapshot has a dura- 
tion of ~1 ks, for a total of up to ~5 ks per day, which is the typical 
time-scale for the flare durations observed in SFXTs (Sguera et al. 
2006). For this reason, we cannot establish the number of observed 
flares and their duration; therefore, we cannot use this information 
to compare the observed light curve with that calculated with the 
clumpy wind model of Ducci et al. (2009), as has been done in 
previous work (Ducci et al. 2009, 2010; Romano et al. 2010). How- 
ever, it is possible to compare the observed cumulative luminosity 
distribution out of eclipse with that calculated with the clumpy wind 
model. For the former, the Swty?/XRT count rate, measured in the 
0.3-10 keV, was converted to the 0. 1-100 keV luminosity using the 
average spectral parameters obtained by Walter et al. (2006) and 
Ducci et al. (2010), and assuming a distance of 13 kpc. For the lat- 
ter, we adopted the physical parameters for the 08.51 primary and 
the neutron star secondary reported above, and the force multiplier 
parameters (see Ducci et al. 2009) obtained by Shimada et al. (1994) 
for an 08.51 star (k = 0.375, a = 0.522, S = 0.099). 

Fig. 6 shows the Swift/XKT cumulative luminosity distribution 
(black solid line) with its 90 per cent uncertainties (black dashed 
lines). Fig. 6 also shows the cumulative luminosity distribution 
(red dot-dashed line) calculated with the clumpy wind model by 
adopting the following parameter values: mass-loss rate M tot = 
6 x 10' 7 MQyr' 1 , terminal velocity Voc = 1000 kms' 1 , /3 = 1, 
fraction of mass lost in clumps / = 0.9, mass distribution power- 
law index f = 1.1, power-law index of the initial clump dimension 
distribution y = — 1, minimum clump mass M a = 5 x 10 18 g and 
maximum clump mass Mi, = l() 20 g. Further acceptable solutions 
can be found assuming the wind parameters in the ranges reported 
in Table 3. 

The wind parameters obtained for IGR J 1 64 18 — 4532 are roughly 
in agreement with those of HMXBs previously studied (see Ducci 
et al. 2009, 2010; Romano et al. 2010), with the exception off 
and Voc ■ The obtained low value of the terminal velocity can be 
explained as follows: when the wind is highly ionized by the X-ray 
photons emitted by the neutron star, the high ionization modifies 
the dynamics of the line-driven stellar wind of OB supergiants 
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Figure 6. Observed cumulative luminosity distribution (solid black line) 
and relative uncertainties (dashed black lines), compared with a possible 
solution (dot-dashed red line; see text). 


Table 3. Wind parameters values for IGR J 16418— 4532. 


Parameter 

Min. value 

Max. value 

Mot 

4 x lO^ 7 Mg yr' 1 

8 x 10- 7 Mq yr -1 

/ 

0.85 

0.95 

^oo 

800 km s^ 1 

1300kms _1 

p 

0.8 

1.2 

$ 

0.8 

1.6 

y 

-1.5 

-0.5 

M a 

5 x 10 16 g 

5 x 10 19 g 

M h 

5 x 10 19 g 

10 21 g 


(especially in close binary systems), leading to a reduction of the 
wind velocity in the direction of the neutron star (see e.g. Stevens & 
Kallman 1990, and references therein). The mass-loss rate is lower 
than the typical mass-loss rate of 08.51 stars, which is of the order 
of «2 x 10~ 6 Mq yr~* (Vink, de Koter & Lamers 2000), but in 
agreement with the hypothesis that the mass-loss rates derived from 
homogeneous-wind model measurements with the Her method are 
overestimated by a factor of 2-10 if the wind is clumpy (see e.g. 
Lepine & Moffat 2008). However, we cannot exclude the possibility 
that the X-ray variability of IGR J 16418 — 4532 is not totally due to 
the accretion of an inhomogeneous wind. In fact, other mechanisms 
could be at work, like transient accretion discs (Taam, Brown & 
Fryxell 1988) and intermittent accretion flow on to a neutron star 
(Lamb et al. 1977; Ducci et al. 2010). 

Our results, based on our intense monitoring of IGR 
J 1 6418 — 4532, show that its X-ray properties (dynamical range, 
flaring activity and peak luminosity) are consistent with an inter- 
mediate SFXT nature. 
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